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The strength distributions of the isoscalar giant monopole resonance have been calculated 
in 16 0, 40 Ca, 90 Zr, 112 ~ 124 Sn, 144 Sm, and 208 Pb nuclei within the self-consistent random phase 
approximation and its extensions which include pairing correlations and quasiparticle-phonon cou- 
pling. The results are compared with the available experimental data. The problem of the nuclear 
matter incompressibility is discussed. 

Theoretical description of the isoscalar giant monopole resonance (ISGMR) is 
an important issue, first of all, because its energy is closely related to the value of the 
incompressibility modulus of infinite nuclear matter (INM) (see [HE]), which 
in turn is a universal characteristic of the effective nuclear forces. The most widely 
used method of evaluating Koo is based on the self-consistent Hartree-Fock (HF) or 
random phase approximation (RPA) calculations of the mean energies of the ISGMR 
using effective Skyrme or Gogny forces (see, e.g., Refs. [H El El IH [5]). Because 
can be calculated from the known parameters of the given force, its value is 
£T) . estimated as the one corresponding to the force that gives the best description of 

the experimental energies. The non-relativistic estimates obtained in such a way 
lead to the value = 210 ± 30 MeV (see, e.g., Refs. [2] E]), though the recent 
results testify to the upper limit of this estimate (see [HE]). 

However, neither the HF approximation nor the RPA can provide full descrip- 
tion of the experimental data. In particular, they cannot reproduce the total width 
T of the giant resonance, since in these models the mechanism responsible for the for- 
mation of the spreading width is absent. At the same time, it is well known that 
the spreading width is a considerable part of T. In the present work we report the 
results of the ISGMR calculations within three models: (i) the RPA; (ii) the quasi- 
particle RPA (QRPA); (iii) the quasiparticle time blocking approximation (QTBA, 
see [HI E]) which is an extension of the QRPA including quasiparticle-phonon cou- 
pling (QPC). The QPC provides the necessary mechanism producing the spreading 
width in the QTBA. The pairing correlations are taken into account both in the 
QRPA and in the QTBA. 
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The basic equation of our approach is the equation for the effective response 
function R (u). It has the same form both in the (Q)RPA and in the QTBA. In 
the shorthand notations it reads (we will follow notations of Ref. [7J) 

R eS {co) = A(u) - A(u) FR^itu) , (1) 

where A{oj) is a correlated propagator and T is an amplitude of the effective resid- 
ual interaction. In the present work we use the version of the QTBA in which 
the ground-state correlations caused by the QPC are neglected. In this case the 
correlated propagator A(u) is defined by the equation 

A{w)=A(u))-A{u))$(u;)A(u), (2) 

where A(uS) is the uncorrelated QRPA propagator, 

= $( rcs )( w ) -$( rcs )(0), (3) 

and <&( res )(u;) is a resonant part of the interaction amplitude responsible for the QPC 
in our model (see Refs. [6j [7J for details). In the case of the (Q)RPA, the amplitude 

(uj) is set equal to zero, and A(u) reduces to the uncorrelated propagator A(u). 

The dynamical pairing effects (particle-particle channel contributions) are taken 
into account both in the QRPA and in the QTBA. This enables one to eliminate 
the + spurious states on the QRPA level. The problem of the + spurious states 
in the QTBA is solved with the help of a combination of the so-called subtraction 
procedure (see Eq. ([3]) and Refs. [6j UJ for details) and the projection technique 
described in [8]. 

The strength function of the ISGMR S(E) is determined by R eS (u) via the 
formulas 

S(E) = — ImU(E + iA), (4) 

7T 

nM = -^Tr((eV°)t J R c V)(eV°)), (5) 

where n(co>) is the nuclear polarizability, E is an excitation energy, A is a smearing 
parameter, V° is an external field, and e is an effective charge operator. In the 
case of the isoscalar + excitations the one-body operator eV° is proportional to the 
identity matrices both in the spin and in the isospin indices. Its radial dependence 
is taken in our calculations in the form eV° = r 2 . The smearing parameter is taken 
to be equal to 500 keV in all the calculations presented here. 

We used a self-consistent calculation scheme based on the HF and Bardeen- 
Cooper-Schrieffer (HF+BCS) approximation. The details of this scheme are de- 
scribed in Ref. [9]. Within the RPA the self-consistency is full. The single-particle 
continuum is also exactly included on the RPA level. The self-consistent mean 
field and the effective interaction (including the spin-orbital and the Coulomb con- 
tributions in both quantities) are derived from the Skyrme energy functional. In 
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the calculations, two Skyrme force parametrizations with the nucleon effective mass 
m* = m were used (see Ref. [ID]): T5 force with = 202 MeV and T6 force with 
Kqo = 236 MeV. The ISGMR calculations were performed in 16 0, 40 Ca, and 90 Zr 
nuclei, in the chain of the Sn isotopes, and in 144 Sm and 208 Pb nuclei. 

The results for the ISGMR strength distributions in the even-A 112_124 Sn iso- 
topes are presented in Fig. Q] and in Table [TJ The mean energies of the ISGMR 
drawn in the tables are defined via the ratios of the energy-weighted moments 
determined as 



The energy interval limited by E\ = 10.5 MeV and E% = 20.5 MeV was taken the 
same as in Ref. [11], where the experimental data on the strength distributions of the 
ISGMR in the tin isotopes were reported. The peak energies E QMR and the widths 
T of the ISGMR were obtained from the Lorentzian fit of the calculated functions 
S(E). As can be seen from the Tabled! the agreement of the theoretical results with 
the experimental mean and peak energies in the case of the T5 Skyrme force is fairly 
good both in the QRPA and in the QTBA. The fact that the mean and peak energies 
obtained in the QRPA and in the QTBA are very close to each other is explained by 
the subtraction procedure used in our calculations (see Eq. (J3j) and Ref. [7] for the 
discussion). The main reason of the agreement with the experiment in this case is 
comparatively low value of the incompressibility modulus of INM (K^ = 202 MeV) 
produced by the T5 Skyrme-force parametrization. The other parametrizations with 

around 240 MeV give too large mean energies of the ISGMR in the considered 
tin isotopes as compared with the experiment. 

For comparison, in Table [1] we draw the QRPA results obtained with the T6 
Skyrme force (i^oo = 236 MeV). As can be seen, the T6 peak energies E QMR are 
greater than the experimental values for the tin isotopes by 1.2-1.9 MeV. This fact 
agrees with the results of Ref. [T2] where it was obtained that the relativistic RPA 
calculations based on the force with = 230 MeV consistently overestimate the 
centroid energies of the ISGMR in the same tin isotopes. 



In Tabled the (Q)RPA and QTBA results for 16 0, 40 Ca, 90 Zr, 144 Sm, and 208 Pb 



nuclei are listed in comparison with the available ISGMR data. As can be seen from 
this table, the T5 Skyrme force leads to a better description of the experiment in 
the light nuclei 16 and 40 Ca. The experimental value of the ISGMR peak energy 
in 90 Zr obtained by the RCNP (Osaka University, Ref. [IT]) and the value mi/m 
for this nucleus obtained by the Gaussian fit in the TAMU experiment (Texas A&M 
University, Ref. [T5]) are also fairly well reproduced by the calculations with the T5 
Skyrme force. The other experimental data for 90 Zr and 144 Sm nuclei lie between 
the values calculated with the T5 and T6 forces. 

It is worth noting that the value = 202 MeV corresponding to the T5 
Skyrme force lies within the interval 210 ± 30 MeV which was considered for a 
long time as the non-relativistic estimate for this quantity. The recent results [H [5] 
inferring to be 230-240 MeV were obtained within the RPA and the constrained 
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Energy (MeV) 

Figure 1: Isoscalar giant monopole resonance in the even- A 112 ~ 124 Sn isotopes cal- 
culated within the QRPA (dashed line) and the QTBA (solid line) using T5 Skyrme 
force. The smearing parameter A is equal to 500 keV. Experimental data (solid 
squares) are taken from Ref. 
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Table 1: Mean energies and Lorentzian-fit parameters for the ISGMR strength dis- 
tributions in the even- A 112_124 Sn isotopes. Theoretical results are obtained within 
the QRPA and the QTBA using T5 and T6 Skyrme forces (indicated in parentheses). 
The mean energies are calculated for 10.5-20.5 MeV energy interval. Experimental 
values are taken from Ref. [IT] (RCNP, Osaka University). 



Method ^/rrii/m-i mi/ mo ^m^jm\ i? GMR T 
(MeV) (MeV) (MeV) (MeV) (MeV) 



Sn QRPA (T6) 17.0 17.1 17.3 17.3 1.9 

QRPA (T5) 15.8 15.9 16.1 15.9 1.8 

QTBA (T5) 15.7 15.8 16.2 15.8 3.7 

Exp. 16.llgj 16.2+g;^ 16.7+°| 16.llgj 4.0±g;^ 

Sn QRPA (T6) 16.9 17.0 17.2 17.3 2.0 

QRPA (T5) 15.7 15.8 16.0 15.8 1.8 

QTBA (T5) 15.6 15.7 16.1 15.7 3.7 

Exp. 15.9±gj 16.llg;l 16.5±g;l 15.9±g;j 4.1±g;l 

Sn QRPA (T6) 16.8 16.9 17.1 17.2 2.1 

QRPA (T5) 15.6 15.6 15.9 15.7 1.9 

QTBA (T5) 15.5 15.6 16.0 15.6 3.8 

Exp. 15.7±gj 15.8±g;j 16.3±g;l 15.8±g;j 

Sn QRPA (T6) 16.6 16.7 17.0 17.1 2.1 

QRPA (T5) 15.4 15.5 15.8 15.6 2.0 

QTBA (T5) 15.4 15.5 15.9 15.5 3.9 

Exp. 15.6+gj 15.8±g;j 16.3±g;j 15.6±gj 4.3±g;t 

Sn QRPA (T6) 16.5 16.6 16.9 17.0 2.2 

QRPA (T5) 15.3 15.4 15.7 15.5 2.1 

QTBA (T5) 15.3 15.4 15.8 15.3 3.9 

Exp. 15.5±g;j 15.7+g;} 16.2±g;l 15.4+g;l 4.9±g;| 

Sn QRPA (T6) 16.4 16.5 16.8 16.9 2.3 

QRPA (T5) 15.2 15.3 15.5 15.4 2.1 

QTBA (T5) 15.1 15.3 15.7 15.2 3.8 

Exp. 15.2±g;j 15.4±g;j 15.9+g;| 15.0+g;l 4.4±g;t 

Sn QRPA (T6) 16.2 16.4 16.7 16.7 2.3 

QRPA (T5) 15.0 15.1 15.4 15.2 2.2 

QTBA (T5) 15.0 15.2 15.5 15.1 3.8 

Exp. 15.118:1 15.3ig;{ 15-8±g;i 14.8+g;l 4.5+g;I 
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Table 2: The same as in Table EEJ but for 16 0, 40 Ca, 90 Zr, 144 Sm, and 208 Pb nuclei. 
The theoretical mean energies are obtained for the following energy intervals: 11-40 
MeV in 16 0, 10-55 MeV in 40 Ca, and 5-25 MeV in 90 Zr, 144 Sm, and 208 Pb. Experi- 
mental data are taken from Refs. [133 [H E2 EE] (TAMU, Texas A&M University) 
and [TTJEB] (RCNP, Osaka University). The data from Ref. [15J for 90 Zr denoted as 
[T5] a and (15] b correspond to the slice analysis and to the Gaussian fit, respectively. 



Method \Jm\jm-\ mi/m ^m^/mi £ GMR T 
(MeV) (MeV) (MeV) (MeV) (MeV) 



3 RPA (T6) 22.7 23.3 25.1 22.0 10.3 

RPA (T5) 21.3 21.8 23.4 20.5 7.2 

Exp. [El 19 1 63g| 2L13g| 24.89±g;jg SJgjj 

Ca RPA (T6) 21.2 21.6 23.0 21.6 4.8 

RPA (T5) 19.7 19.9 21.1 19.7 4.0 

Exp. M 19.18lgf 4.88±g;i; 

'Zr QRPA (T6) 18.0 18.2 18.6 18.0 3.0 

QTBA (T6) 18.0 18.2 18.7 18.0 3.8 

QRPA (T5) 16.6 16.8 17.2 16.5 2.0 

QTBA (T5) 16.6 16.8 17.3 16.5 2.6 

Exp. [15] a 17.8lt°;| 17.89i°;|° 

Exp. US] b 16.80 

Exp. gH 16.6lgj 4.9l°; 2 

Sin QRPA (T6) 15.8 16.0 16.4 15.8 2.0 

QTBA (T6) 15.8 16.1 16.7 15.8 3.0 

QRPA (T5) 14.6 14.7 15.1 14.5 1.5 

QTBA (T5) 14.6 14.8 15.4 14.4 2.5 

Exp. m 15.40i°2j 3.40l8:ig 

Exp, m i5.3oi°:l2 3.7118:63 

*Pb RPA (T6) 13.8 14.0 14.5 13.9 1.9 

QTBA (T6) 13.8 14.1 14.8 13.9 3.2 

RPA (T5) 12.6 12.7 13.2 12.6 1.6 

QTBA (T5) 12.6 12.8 13.5 12.5 2.8 

Exp. US] 13.9618J8 2.881818 

Exp. m 13-4181 4.0181 



HF method on the base of the experimental data in fact only for the one nucleus 
208 Pb. The mean energy mi/m of the ISGMR in 208 Pb obtained in the TAMU 
experiment [16J is also nicely reproduced in our calculations when the T6 force with 
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Koo = 236 MeV is used. On the other hand, the T5 force gives the QTBA value 
mi/mo for 208 Pb which is lesser by 1.2 MeV as compared with this experiment. 
Note, however, that the ISGMR data even for the well-studied nucleus 208 Pb are 
not so far quite unambiguous. In particular, the experimental value of the ISGMR 
peak energy in 208 Pb measured in the RCNP experiment [17] is lesser by 0.6 MeV as 
compared with the value mi /mo from Ref. [16] and lies between the QTBA values 
of E QMR obtained with the T5 and T6 Skyrme forces. 

In contrast to the mean energies, the (Q)RPA and the QTBA give substantially 
different values for the total width T of the ISGMR. As was mentioned above, the 
reason is that the (Q)RPA does not produce the spreading width T^, whereas in 
the QTBA it is formed due to the QPC effects. The fact that is a considerable 
part of T is illustrated by the results listed in the tables: the (Q)RPA strongly 
underestimates the experimental values of the total width, while reasonably good 
agreement is achieved in the QTBA. The exception is the case of the light nuclei in 
which the resonance width is fairly well reproduced already within the RPA. 

In conclusion, we have presented the results of the theoretical analysis of the 
ISGMR strength distributions in 16 0, 40 Ca, 90 Zr, 112 ~ 124 Sn, 144 Sm, and 208 Pb nu- 
clei. The calculations were performed within the (Q)RPA and the QTBA which 
is an extension of the QRPA including quasiparticle-phonon coupling (QPC). We 
used calculational scheme based on the HF+BCS approximation which is fully self- 
consistent on the RPA level. In the calculations, two Skyrme force parametrizations 
were used. The T5 parametrization with comparatively low value of the incompress- 
ibility modulus of infinite nuclear matter (K^ = 202 MeV) allowed us to achieve 
good agreement with the experimental data for tin isotopes within the QTBA in- 
cluding resonance widths. However, this parametrization underestimates the ex- 
perimental ISGMR mean energy for the 208 Pb nucleus which is usually used in the 
fit of the Skyrme force parameters. On the other hand, the T6 Skyrme force with 

= 236 MeV nicely reproduces the ISGMR mean energy for 208 Pb but overesti- 
mates the energies for 112_124 Sn isotopes by more than one MeV. For the light nuclei 
16 and 40 Ca, the T5 Skyrme force gives a better description of the experimen- 
tal ISGMR mean energies as compared with the T6 force. The experimental data 
on the ISGMR energies in 90 Zr and 144 Sm nuclei lie between the values calculated 
by us with the T5 and T6 forces. On the whole, these results do not allow us to 
decrease the ambiguity in the value of Koo as compared with the known estimate 
Kqo = 210 ±30 MeV. However, our calculations confirm that an account of the QPC 
is necessary to describe the total width of the ISGMR. 
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